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Abstract 

' In a combined analysis of the experimental data on the coupled processes vrvr 

. TTTT^KK in the channel with I^J^'-^ = O'^O'*"^, the following results are obtained in 

a model-independent way: 1) The /o(665) state with properties of the c-meson is 
proved to exist; 2) it is shown that the /o(980) and especially /o(1370) (if exists) have 
. a dominant ss component; 3) indications of the glueball nature of the /o(1500) and of 

a considerable ss component in the /o(1710) are obtained; 4) conclusion on the linear 
■ realization of chiral symmetry is drawn. 

^ — I 

^ ■ 1 Introduction 

O I The problem of scalar mesons is far off the solution up to now. For instance, at present, 

Q I additional arguments have been added by N.N. Achasov in favour of the 4-quark nature 

CN ■ of /o(980) and ao(980) mesons on the basis of interpretation of the experimental data on the 

^ decays (j) 777^71°, "yTT^rj 0. On the other hand, F.E. Close and A. Kirk have shown that 



O 



mixing between the /o(980) and ao(980) radically affects some existing predictions of their 
produciton in radiative decay. Generally, the 4-quark interpretation, beautifully solving 
^1 the old problem of the unusual properties of scalar mesons, sets new questions. Where are 
i-^H ■ the 2-quark states, their radial excitations and the other members of 4-quark multiplets 

9, 9*, 36 and 36*, which are predicted to exist below 2.5 GeV @]? 
^ ! The discovered states in the scalar sector and their properties do not allow one to make up 

^ I the scalar qq nonet and to solve other exiting questions up to now. To draw this conclusion, 
it is sufficient to compare various variants of assignment of these states to quark-model con- 
figurations, for example, Generally, difficulites in understanding the scalar-isoscalar 



sector seem to be related to both the hard-accounting influence of the vacuum (and such 
effects as the instanton contributions) and a strong model dependence of an information 
about wide multichannel states. 

Obviously, it is important to have a model-independent information on investigated states 
and on their QCD nature. It can be obtained only on the basis of the first principles 
(analyticity, unitarity) immediately applied to experimental data analysis. Earlier, we have 
proposed this method for 2- and 3-channel resonances and developed the concept of standard 
clusters (poles on the Riemann surface) as a qualitative characteristic of a state and a 
sufficient condition of its existence IT^. We outline this below for the 2-channel case of the 
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coupled processes tttt tttt, KK in the channel with J^*-^ = 0+0++. Since, in this work, 
a main stress is laid on studying lowest states, it is sufficient to restrict oneself to a two- 
channel approach when considering simultaneously the coupled processes tttt tttt, KK, 
though in the future it is nesessary to take into account the thresholds of other coupled 
processes, ffist of all, of rjrj and rjr]' scatterings. In this work, we are going to show that 
the large background, which one has obtained earlier in various analyses of the s-wave tttt 
scattering [T^, hides, in reality, the a- meson []T^ below 1 GeV and the effect of the left-hand 
branch-point. Therefore, in the uniformizing variable, one must take into account, besides 
the branch-points corresponding to the thresholds of the processes tttt — > tttt, KK, also the 
left-hand branch-point at s = 0, related to the background in which the crossing-channel 
contributions are contained. Furthermore, we shall obtain definite indications of the QCD 
nature of other /q resonances and of the linear realization of chiral symmetry. In part the 
results of this report have been published in the works ||T5[| . 

Note that recent new analyses of old and new experimental data found a candidate 
for the cr-meson below 1 GeV (see, e.g., 0, |T^-|2^). However, these analyses use 
either the Breit - Wigner form that is insufficiently- flexible even if modified or the /T- matrix 
formalism without taking account of the energetic-closed (maybe, important) channels, or 
specific forms of interactions in the quark models; therefore, there one cannot talk about 
the model independence of results. Besides, in these analyses, a large vrvr-background is 
obtained. 

The layout of this report is as follows. In Section II, we outline the two-coupled-channel 
formalism, determine the pole clusters on the Riemann surface as characteristics of multi- 
channel states, and introduce a new uniformizing variable, allowing for the branch-points 
of the right-hand (unitary) and left-hand cuts of the Trvr-scattering amplitude. In Section 
III, we analyze simultaneously experimental data on the processes tttt — > tttt, KK in the 
isoscalar s-wave on the basis of the presented approach. Note that earlier we have shown in 
two approaches and ||T5[ without and with taking into account the left-hand branch-point 
^/s in the uniformizing variable, respectively, that the minimal scenario of the simultaneous 
description of two coupled processes vrvr tth, KK is realized without the /o(1370) and 
/o(1710) that are in the Particle Data Group tables. Therefore, we consider also the variants 
of description of the indicated processes including these states separately as well as simul- 
taneously, and obtain indications of their QCD nature different from the results of many 
other works (note that our approach is based on the first principles and is free from dynamic 
assumtions, therefore, our results are rather model-independent). In the Conclusion, the 
obtained results are discussed. 



2 Two-Coupled-Channel Formalism 

Here we restrict ourselves to a 2-channel consideration of the coupled processes tttt 
TTTT, KK. Therefore, we have the 2-channel ^-matrix determined on the 4-sheeted Riemann 
surface. The matrix elements Saf^, where a,j3 = l{7r7r),2{KK), have the right-hand cuts 
along the real axis of the s-plane, starting at 4m^ and ^m\, and the left-hand cuts, begin- 
ning at s = for 5*11 and at 4(m|^ — m^) for S22 and 512. The Riemann-surface sheets are 
numbered according to the signs of analytic continuations of the channel momenta 

h = {s/A - mlfl\ k2 = {s/A - mlfl' (1) 

as follows: signs (Im/ci, ImA;2) = ++, — h, , H — correspond to the sheets I, II, HI, IV. 
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To obtain the resonance representation on the Riemann surface, we express analytic 
continuations of the matrix elements to the unphysical sheets S"^^ {L = II, III, IV) in 
terms of those on the physical sheet S'^^: 
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Here det = ^(^^2^2 - ('^(2)^ ('^(2)^ = s'^^is - Aml){s - 4ml)F{s); in the limited 
energy interval, F{s) is proportional to the squared product of the coupling constants of the 
considered state with channels 1 and 2. These formulas are convenient by that the S'-matrix 
elements on the physical sheet S^j^ have, except for the real axis, only zeros corresponding 
to resonances, at least, around the physical region that is interesting for us. Formulas (H) 
immediately give the resonance representation by poles and zeros on the 4-sheeted Riemann 
surface. One must distinguish between three types of 2-channel resonances described by a 
pair of conjugate zeros on sheet I: (a) in Su, (b) in 6*22, (c) in each of and 822- As seen 
from eqs. (0), to the resonances of types (a) and (b), there corresponds a pair of complex 
conjugate poles on sheet III shifted relative to a pair of poles on sheet II and IV, respectively. 
For the states of type (c), one must consider the corresponding two pairs of conjugate poles on 
sheet III. A resonance of every type is represented by a pair of complex-conjugate clusters (of 
poles and zeros on the Riemann surface) of a size typical of strong interactions. The cluster 
kind is related to the state nature. The resonance coupled relatively more strongly to the vrvr 
channel than to the KK one is described by the cluster of type (a); in the opposite case, it is 
represented by the cluster of type (b) (say, the state with the dominant ss component); the 
flavour singlet {e.g. glueball) must be represented by the cluster of type (c) as a necessary 
condition, if this state lies above the thresholds of considered channels. 

Furthermore, according to the type of pole clusters, we can distinguish, in a model- 
independent way, a bound state of colourless particles {e.g., KK molecule) and a qq bound 



state |12, 21]. Just as in the 1-channel case, the existence of a particle bound-state means 
the presence of a pole on the real axis under the threshold on the physical sheet, so in the 
2-channel case, the existence of a particle bound-state in channel 2 {KK molecule) that, 
however, can decay into channel 1 (tttt decay), would imply the presence of a pair of complex 
conjugate poles on sheet II under the second-channel threshold without an accompaniment 
of the corresponding shifted pair of poles on sheet III. Namely, according to this test, earlier, 
the interpretation of the /o(980) state as a KK molecule has been rejected. 

For the simultaneous analysis of experimental data on coupled processes, it is convenient 
to use the Le Couteur- Newton relations |2^ expressing the ^-matrix elements of all coupled 
processes in terms of the Jost matrix determinant d{ki, ^2), the real analytic function with 
the only square-root branch-points at fcj = 0. To take into account, in addition to the latter. 
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also the left-hand branch-point at s = 0, the uniformizing variable is used [| 



(3) 



It maps the 4-sheeted Riemann surface with two unitary cuts and the left-hand cut onto the 
v-plane. (Note that other authors have used the parameterizations with the Jost functions at 



analyzing the s-wave tttt scattering in the one-channel approach |Q and in the two-channel 
one [^. In latter work, the uniformizing variable /c2 has been used, therefore, their approach 
cannot be employed near by the tttt threshold.) In Fig.|l], the plane of the uniformizing 









Im V 




I 








I 


© 






i 








is- 








-1 1 




1 b 


Rev 






<IV 






in 








in 



Figure 1: Uniformization plane for the vrvr-scattering amplitude. 

variable v for the Trvr-scattering amplitude is depicted. The Roman numerals (I,. . . , IV) 
denote the images of the corresponding sheets; the thick line represents the physical region; 
the points i, 1 and b = {rriK + m,^) / {jtik — rn.,,) correspond to the tttt, KK thresholds and 
s = oo, respectively; the shaded intervals (— oo, — 6], [— [6, cxd) are the images of 
the corresponding edges of the left-hand cut. The depicted positions of poles (*) and of zeros 
(o) give the representation of the type (a) resonance in 5*11. 

On the f-plane, Sn has no cuts; however, 5*12 and 6*22 do have the cuts which arise from 
the left-hand cut on the s-plane, starting at s = 4:{m\ — m'^), which further is neglected in 
the Riemann- surface structure, and the contribution of this cut is taken into account in the 
KK background as a pole on the real s-axis on the physical sheet in the sub-ii'-ft'-threshold 
region. 



On f-plane, the Le Couteur-Newton relations are [|T2|, ^2 



S 
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s. 



d{v- 



d{v) ' d{v) ' '^'"^'^ d[v) ■ 

The (i(t>)-function already does not possess branch-points and is taken as 

d = dsdres, (5) 

^The analogous uniformizing variable has been used, e.g., in Ref. [2^ at studying the forward elastic pp 
scattering amplitude. 



SiiS- 
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d{-v) 



(4) 
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where = B^^Bk] B^^ contains the possible remaining Trvr-background contribution, related 
to exchanges in crossing channels (the consequent analysis gives 5^ = 1); Bx is that part of 
the KK background which does not contribute to the Trvr-scattering amplitude: 

Bk = v-\1-vov)\1 + v*v)\ (6) 

The fourth power in is explained in our work [|I^. The function dresi^) represents the 
contribution of resonances, described by one of three types of the pole-zero clusters, i.e., 

M 

dres = V-'''l[il-v:v)il + VnV), (7) 
n=l 

where M is the number of pairs of the conjugate zeros. 



3 Analysis of experimental data 



We simultaneously analyze the available experimental data on the vrvr-scattering [^] and 
the process tttt — > KK H in the channel with I^J^^ = 0+0++. As the data, we use the 
results of phase analyses which are given for phase shifts of the amplitudes (5i and ^12) and 
for moduli of the S'-matrix elements t]i = \Saa\ {a.=l'KTT ,2K K) and ^ = |5'i2|. The 2-channel 
unitarity condition gives i]i = r]2 = f]^ ^ = (1 ~ v'^Y^'^j ^12 = 5i + 62- 

We consider four variants, in which the following states are taken into account: 
Variant 1 : The /o(665) and /o(980)) with the clusters of type (a), and /o(1500), of type (c); 
The same three resonances + the /o(1370) of type (b); 
The /o(665), /o(980)) and /o(1500) + the /o(1710) of type (b); 
All the five resonances of the indicated types. 



Variant 2 



Variant 3 



Variant 4 




Figure 2: 
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Table 1: 



Variant 


1 


2 


Quantity 




5i2 e 






Number of 
exp. points 


160 50 


80 33 


160 50 


80 39 


2 /TVT 


2.7 0.72 


2.37 1.1 


2.85 0.82 


3.98 0.92 




1.96 


2.01 


2.01 


3 


2 /AT 


1.98 


2.45 


(so, GeV ) 


0.954381 + 0.29859i 
(0.441) 


0.97925 + 0.202657i 
(0.6466) 


Variant 


3 


4 


Quantity 




5i2 e 






Number of 
exp. points 


160 50 


80 42 


160 50 


80 42 


xVNdf 


2.38 0.8 


2.25 0.92 


2.57 0.85 


4.74 1.05 


xVNdf 


1.72 


1.8 


1.81 


3.49 




1.76 


2.59 


(so,GeV2) 


0.954572 + 0.29798i 
(0.4646) 


0.982091 + 0.188405i 
(0.678) 



The other possibihties of the representation of these states are rejected by our analysis. 
We consider these variants, because the minimal possibility of the simultaneous description 
of two coupled processes tttt 'kh,KK is realized without the /o(1370) and /o(1710) [|15[] . 

The TTTT-scattering data are described from the threshold to 1.89 GeV in all the four 
variants (in addition, we take B^^ = 1) and are taken from the analysis by B. Hyams et 
al. in this energy region, and below 1 GeV, from many works [^. For the reaction 
TTTT — s> KK^ practically all the accessible data are used, but the description ranges are slightly 
different for various variants and extend from the threshold to ~ 1.4 GeV for variant 1, to 
~ 1.46 GeV for variant 2, and to ~ 1.5 GeV for variants 3 and 4. Table |l| demonstrates the 
quality of fits to the experimental data (the number of fitted parameters is 17 for variant 1, 
21 for variants 2 and 3, 25 for variant 4). When calculating x^/Ndf, we have rejected the 
experimental points at 0.61, 0.65, and 0.73 GeV for 5i, at 0.99, 1.65, and 1.85 GeV for rj, at 
1.111, 1.163, and 1.387 GeV for 5i2, and at 1.002, 1.265, and 1.287 GeV for ^ that give an 
especially large contribution to y^. We note that two variants (1 and 3) are the best, both 
without the /o(1370), and we stress that this analysis uses the parameterless description of 
the Tin background. 

Let us indicate the obtained zero positions, on the v plane, of the corresponding reso- 
nances best variant 3: 



for/o(665) : 


Vl 


= 1.38633 + 0.230588i, 


V2 = 


0.904085 - 


- 0.263033i, 


for/o(980) : 


V3 


= 1.05103 + 0.0487473Z, 


= 


0.864109 - 


- 0.0922272^, 


for/o(1500) : 


V5 


= 1.2477 + 0.0321349i. 


V6 = 


1.24027- 


0.0384191i, 




V7 


= 0.804333 - 0.0179899i, 


Vs = 


0.795579 - 


- 0.02539852, 


for/o(1710) : 


V9 


= 1.25928 - 0.0115127Z, 


VlO = 


= 0.795429 


- 0.00629969Z 
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Table 2: Pole clusters for considered resonances in variant 3. 



Sheet 


II 


III 


IV 


/o(665) 


E, MeV 
r, MeV 


570±13 
590±24 


700±15 
72±5 




/o(980) 


E, MeV 
r, MeV 


989±5 
29±7 


982±14 
195±21 




/o(1500) 


E, MeV 
r, MeV 


1505±23 
272±25 


1490±30 1510±25 
220±26 370±30 


1430±20 
275±32 


/o(1710) 


E, MeV 
r, MeV 




1680±20 
114±15 


1700±15 
138±21 


Table 3: Pole clusters for considered resonances in variant 4. 


Sheet 


II 


III 


IV 


/o(665) 


E, MeV 
r, MeV 


600±16 
605±28 


715±17 
59±6 




/o(980) 


E, MeV 
r, MeV 


985±5 
27±8 


984±18 
210±22 




/o(1370) 


E, MeV 
r, MeV 




1310±22 
410±29 


1320±20 
275±25 


/o(1500) 


E, MeV 
r, MeV 


1528±22 
385±25 


1490±30 1510±20 
220±24 370±30 


1510±21 
308±30 


/o(irio) 


E, MeV 
r, MeV 




1700±25 
86±16 


1700±20 
115±20 



Figures ^ demonstrate the comparison of the obtained energy dependence of four analyzed 
quantities with the experimental data: The short-dashed lines correspond to variant 1; the 
long-dashed curves, to variant 2; the dot-dashed ones, to variant 4; and the solid lines, to 
variant 3. 



Table 4: 





/o(665) 


/o(980) 


/o(1370) 


/o(1500) 


9i, GeV 


0.652 ±0.065 


0.167 ±0.05 


0.116 ±0.03 


0.657 ±0.113 


92, GeV 


0.724 ±0.1 


0.445 ± 0.031 


0.99 ±0.05 


0.666 ±0.15 



In Tables ^ and |], the obtained pole clusters of considered resonances are shown on the 
corresponding sheets on the complex energy plane {^/s^ = E^ — iTr) for the best variant 3 
(without the /(1370)) and for variant 4 (with all five states). 

The coupling constants of obtained states with vrvr (gi) and KK ((72) systems are calcu- 
lated through the residues of amplitudes at the pole on sheet II - for resonances of types 
(a) and (c), and on sheet IV - for resonances of type (b). Expressing the T-matrix via the 
^'-matrix as 

Sii = l + 2ipiTii, S12 = 2iy/pip2Ti2, (8) 
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where pi = ^ (s — Arnf)/ and taking the resonance part of the amphtude as 

= T.9^rgr,D;\s) (9) 

r 

with Dr{s) being an inverse propagator (-Dr(s) oc s — Sr), we show the results of that 
calculation in Table ||. We see that the /o(980) and especially the /o(1370) are coupled 
essentially more strongly to the KK system than to the vrvr one, which tells about the 
dominant ss component in these states. The /o(1500) has the approximately equal coupling 
constants with the vrvr and KK systems, which apparently could point up to its dominant 



glueball component The coupling constant of the /o(1710) with the vrvr channel cannot 
be calculated by this method, unless the description of vrvr KK reaction is obtained in 
the region of this resonance. But this state is represented by the cluster corresponding to 
the dominant ss component. 

Let us also present the calculated scattering lengths: For the KK scattering: 
a° = -1.25 ±0.11 + (0.65 ±0.09)i, [m;|]; (variantl), 
= -1.548 ±0.13 ± (0.634 ±0.1)i, [m;^]; (variant2), 
a° = -1.19 ±0.08 ± (0.622 ±0.07)i, [m;|]; (variant3), 
= -1.58 ±0.12 ± (0.59 ±0.1)i, [m;^]; (variant4). 
Variants 2 and 4 include the /o(1370) unlike variants 1 and 3. We see that Re a^{KK) is 
very sensitive to whether this state exists or not. 

In Table ^ we compare our results for the tttt scattering length with results of some 
other theoretical and experimental works. 



Table 5: 





References 


Remarks 


0.27 ±0.06 (1) 
0.267 ±0.07 (2) 
0.28 ±0.05 (3) 
0.27 ±0.08 (4) 


our paper 


model-independent approach 


0.26 ±0.05 


L. Rosselet et al.[||] 


analysis of the decay K Trvrez/ 
using Roy's model 


0.24 ±0.09 


A. A. Bel'kov et al.pSf 


analysis of 7r~p n^Ti^n 
using the effective range formula 


0.23 


S. Ishida et al.|18 


modified analysis of tttt scattering 
using Breit-Wigner forms 


0.16 


S. Weinberg |8| 


current algebra (non-linear a-model) 


0.20 


J. Gasser, H. Leutwyler |29| 


one-loop corrections, non-linear 
realization of chiral symmetry 


0.217 


J. Bijnens et al.|30(| 


two-loop corrections, non- linear 
realization of chiral symmetry 


0.26 


M.K. Volkov 


linear realization of chiral symmetry 


0.28 


A.N.Ivanov, N.Troitskaya jgl 


a variant of chiral theory with 
linear realization of chiral symmetry 



At first, let us remark about the result of Ref. p8[] for the tttt scattering length. We 
think that such a small value (0.23 m~l) has been obtained, because there has been used an 
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assumption about the negative vrvr background in the phase shift. Now, from Table ^ we see 
that our results correspond to the linear realization of chiral symmetry. 

We have here presented model-independent results: the pole positions, coupling constants 
and scattering lengths. Masses and widths of these states that should be calculated from the 
obtained pole positions and coupling constants are highly model-dependent. For instance, if 
we suppose that the /o(665) is the a-meson, then from the known relation 

gaTTTT = (ml - ml)/V2f^o 

(here f-j^o is the constant of weak decay of the 7r°: /^o = 93.1 MeV), we obtain rrif^ ^ 342 
MeV. 

If we take the resonance part of amplitude in a non-relativistic form ([l^, p. 214) 

j^res _ 

~ E^-E-tT tot/2' 

then we have (variant 3) E„ ^ 570 ± 13 MeV and Ttot ~ 1180 ± 48 MeV; for the so-called 
relativistic form 

rpres _ \/sTel 

ml - S-iy/sTtot' 

the following values are obtained: m„ ^ 820 ± 20 MeV and T ^ 1180 ± 48 MeV. 



4 Conclusions 

On the basis of a simultaneous description of the isoscalar s-wave channel of the pro- 
cesses TTTT TTTT, KK wlth a parameterless representation of the vrvr background, a model- 
independent confirmation of the a-meson below 1 GeV is obtained. We emphasize that 
this is a real evidence of this state, because we have not been enforced to construct the vrvr 
background. 

A parameterless description of the vrvr background is given only by allowance for the 
left-hand branch-point in the proper uniformizing variable. This seems to be related to the 
fact that the exchanges by nearest a- and p-mesons in the crossing channels contribute to 
the TTTT-scattering amplitude with opposite signs and compensate each other. 

Since all the fitted parameters in describing the vrvr scattering are only the positions of 
poles corresponding to resonances, we conclude that our model-independent approach is a 
valuable tool for studying the realization schemes of chiral symmetry. The existence of the 
low-lying state /o(665) with the properties of the a-meson and the obtained vrTr-scattering 
length (aQ(7r7r) f» 0.27[m~+]) suggest the linear realization of chiral symmetry. 

The analysis of the used experimental data gives the evidence that the /o(980) and es- 
pecially /o(1370) resonance (if exists - variants 2 and 4), have the dominant ss component. 
Note that a minimum scenario of the simultaneous description of processes tttt tttt, KK 
goes without the /o(1370) resonance. The best total x^/Ndf for both the analyzed pro- 
cesses is obtained with the set of states: /o(665), /o(980), /o(1500) and /o(1710). The KK 
scattering length is very sensitive to whether the /o(1370) state exists or not. 

The /o(1500) has the approximately equal coupling constants with the tttt and KK 
systems, which apparently could point up to its dominant flavour-singlet {e.g., glueball) 
component . 
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The /o(1710) is represented by the cluster corresponding to the state with the dominant 
ss component. This is in accord with Refs.|[Tl|, ^ where the /o(1500) has been considered 
as a candidate for the scalar glueball. Note that QCD sum rules and the K-matrix 
method showed both the /o(1500) and /o(1710) are mixed states with large admixture 
of the glueball component. Their conclusions about the glueball component concord to our 
conclusion as to the /o(1500) but not /o(1710). It seems that the complement of the combined 
consideration by the rji] and rir]' channels should not change substantially this situation in 
view of the relative distance of the rjr] threshold, and because the glueball component is 



not coupled to the t]T]' system. Note also that the conclusion of QCD sum rules |Q about 
the existence of light glueballs (below 1 GeV) contradicts the lattice calculations and is not 
confirmed by our method. 

We stress that our results are very decisive, because our approach is based only on the first 
principles (analyticity-microcausality and unitarity), immediately applied to the analysis of 
experimental data, and it is free from dynamical assumptions, because a way of its realization 
is based on the mathematical fact that a local behaviour of analytic functions determined 
on the Riemann surface is governed by the nearest singularities on all corresponding sheets. 
It is very important that we were able to describe the considered coupled processes without 
diminishing the number of fitted parameters by some dynamical assumptions. 

We think that multichannel states are most adequately represented by clusters, i.e., by 
the poles on all the corresponding sheets. Pole clusters give a main effect of resonances, 
and on the uniformization plane they are their good representation. The pole positions 
are rather stable characteristics for various models, whereas masses and widths are very 
model-dependent for wide resonances. 

References 

[1] N.N. Achasov, Nucl. Phys. A 675, 279c (2000). 

[2] M.N. Achasov et al, Phys. Lett. B 438, 441 (1998); ibid. 440, 442 (1998). 

[3] F.E. Close, A. Kirk, Phys. Lett. B 515, 13 (2001). 

[4] R.L. Jaffe, Phys. Rev. D 15, 267, 281 (1977). 

[5] N.A. Tdrnqvist, M. Roos, Phys. Rev. Lett. 76, 1575 (1996). 

[6] V.V. Anisovich, D.V. Bugg, A.V. Sarantsev, Phys. Rev. D 58, 111503 (1998). 

[7] W. Lee, D. Weingarten, Phys. Rev. D 59, 094508 (1999). 

[8] D. Black, A.H. Fariborz, J. Schechter, Phys. Rev. D 61, 074001 (2000). 

[9] F.E. Close, A. Kirk, Phys. Lett. B 483, 345 (2000). 

[10] CM. Shakin, Huangsheng Wang, Phys. Rev. D 63, 014019 (2000). 

[11] M.K. Volkov, V.L. Yudichev, Eur. Phys. J. A 10, 223 (2001). 

[12] D. Krupa, V.A. Meshcheryakov, Yu.S. Surovtsev, Nuovo Cim. A 109, 281 (1996); Yad. 
Fiz. 43, 231 (1986); Czech. J. Phys. B 38, 1129 (1988). 



10 



[13] Review of Particle Physics, Europ. Phys. J. C 15, 1 (2000). 

[14] Y. Nambu, G. Jona-Lasinio, Pliys. Rev. 122, 345 (1961); M.K. Volkov, Ann. Phys. 157, 
282 (1984); T. Hatsuda, T. Kunihiro, Phys. Rep. 247, 223 (1994); R. Delbourgo, M.D. 
Scadron, Mod. Phys. Lett. A 10, 251 (1995). 

[15] Yu.S. Surovtsev, D. Krupa, M.Nagy Phys. Rev. D 63, 054024 (2001); Acta Physica 
Polonica B 31, 2697 (2000); Acta Physica Slovaca 51, 189 (2001). 

[16] B.S. Zou, D.V. Bugg, Phys. Rev. D 48, R3948 (1993); ibid. 50, 591 (1994). 

[17] M. Svec, Phys. Rev. D 53, 2343 (1996). 

[18] S. Ishida et al, Progr. Theor. Phys. 95, 745 (1996); ibid. 98, 621 (1997). 
[19] R. Kamihski, L. Lesniak, B. Loiseau, Eur. Phys. J. C 9, 141 (1999). 
[20] L. Li, B.-S. Zou, G.-he Li, Phys. Rev. D 63, 074003 (2001). 
[21] D. Morgan, M.R. Pennington, Phys. Rev. D 48, 1185 (1993). 

[22] K.J. Le Couteur, Proc. Roy. Soc. A 256, 115 (1960); R.G. Newton, J. Math. Phys. 2, 
188 (1961); M. Kato, Ann. Phys. 31, 130 (1965). 

[23] B.V. Bykovsky, V.A. Meshcheryakov, D.V. Meshcheryakov, Yad. Fiz. 53, 257 (1990). 

[24] J. Bohacik, H. Kiihnelt, Phys. Rev. D 21, 1342 (1980). 

[25] B. Hyams et al, Nucl. Phys. B 64, 134 (1973); ibid. 100, 205 (1975); A. Zylbersztejn 
et al, Phys. Lett. B 38, 457 (1972); P. Sondercggcr, P. Bonamy, in Proc. 5th Intern. 
Conf. on Elementary Particles (Lund, 1969), paper 372; J.R. Bcnsinger et al, Phys. 
Lett. B 36, 134 (1971); J.P. Baton et al, Phys. Lett. B 33, 525, 528 (1970); P. Baillon 
et al, Phys. Lett. B 38, 555 (1972); L. Rosselet et al, Phys. Rev. D 15, 574 (1977); 
A.A. Kartamyshev et al, Pis'ma v Zh. Eksp. Tcor. Fiz. 25, 68 (1977); A.A. Bel'kov et 
al, Pis'ma v Zh. Eksp. Teor. Fiz. 29, 652 (1979). 

[26] W. Wetzel et al, Nucl. Phys. B 115, 208 (1976); V.A. Polychronakos et al, Phys. Rev. 

D 19, 1317 (1979); P. Estabrooks, Phys. Rev. D 19, 2678 (1979); D. Cohen et al, Phys. 
Rev. D 22, 2595 (1980); G. Costa et al, Nucl. Phys. B 175, 402 (1980); A. Etkin et 
al, Phys. Rev. D 25, 1786 (1982). 

[27] C. Amsler, F.E. Close, Phys. Rev. D 53, 295 (1996). 

[28] S. Weinberg, Phys. Rev. Lett. 17, 616 (1966). 

[29] J. Gasser, H. Leutwyler, Ann. Phys. 158, 142 (1984). 

[30] J. Bijnens et al, Phys. Lett. B 374, 210 (1996). 

[31] M.K. Volkov, Phys. Elem. Part. At. Nucl. 17, part 3, 433 (1986). 

[32] A.N. Ivanov, N.I. Troitskaya, Nuovo Cim. A 108, 555 (1995). 

[33] S. Narison, Nucl. Phys. B 509, 312 (1998). 



11 



